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Abstract 

The effect of viscosity on the wave resistance experienced by a 
two-dimensional perturbation moving at uniform velocity over the free 
surface of a fluid is investigated. The analysis is based on Rayleigh's 
linearized theory of capillary- gravity waves. It is shown in particular 
that the wave resistance remains bounded as the velocity of the per- 
turbation approaches the minimum phase speed c™" = 
(/9 is the liquid density, 7 is the liquid-air surface tension, and g the 
acceleration due to gravity), unlike what is predicted by the inviscid 
theory. 
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Consider a body of liquid in equilibrium in a gravitational field and having 
a planar free surface. If, under the action of some external perturbation, the 
surface is moved from its equilibrium position at some point, motion will 
occur in the liquid. This motion will be propagated over the whole surface in 
the form of waves, which are called capillary- gravity waves [Q. These waves 
are driven by a balance between the liquid inertia and its tendency, under 
the action of gravity and under surface tension forces, to return to a state 
of stable equilibrium. For an inviscid liquid of infinite depth, the relation 
between the circular frequency u and the wave number k (i.e., the dispersion 
relation) is given by u"^ = gk + 7^^/ p where p is the liquid density, 7 the 
liquid-air surface tension, and g the acceleration due to gravity 0. The 
above equation may also be written as a dependence of wave speed c = u/k 
on wave number 

c=ig/k + ^k/p)'/' (1) 

An important feature of Eq. (|l]) is that it implies a minimum phase speed 
of c™" = {Ag^/pyl'^ reached at fcmm = ^ where = [l/{pg)f^'^ is the 
capillary length For water with 7 = 73 mN m~^ and p = 10^ kg m~^, the 
minimun phase speed is c™" = 0.23 m s~^ and the corresponding wavelength 
is A™" = 271 /k = 1.7 10~^ m. The dispersive property of capillary-gravity 
waves is responsible for the complicated wave pattern generated at the free 
surface of a still liquid by a disturbance moving with a velocity V greater 
than c™™ [0. The disturbance may be produced by a small object (such as a 
fishing line) immersed in the liquid or by the application of an external surface 
pressure distribution Pext- The waves generated by the moving disturbance 
continuously remove energy to infinity. Consequently, for V > c*"*", the 
disturbance will experience a drag, R, called the wave resistance For 
V < c"^*", the wave resistance is equal to zero since, in this case, no waves 
are generated by the disturbance. A few years ago 0, it has been predicted 
that the wave resistance corresponding to a surface pressure distribution 
symmetrical about a point should be discontinuous at = c"^*" [|]. This 
prediction has been checked very recently by Browaeys and co-workers using 
a magnetic fluid 0. The experimental results of Browaeys et al. indicate, 
however, that the disturbance experienced a small but nonzero drag for V < 
c™*". Since this nonzero drag might be due, in part, to the finite viscosity 
of the fluid, it is of some importance to incorporate this physical parameter 
in the inviscid model of Ref. In order to simplify the discussion, we will 
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here consider the case of a pressure distribution Pext locahzed along a hne. 
The more comphcated case of an axisymmetric surface pressure distribution 
will be consider elsewhere. 

We take the xy-plane as the equilibrium surface of the fluid and assume 
that a pressure distribution of the form 

Pext = Pq I 2\ 

travels over the surface in the x-direction with a velocity V (in all what 
follows we assume that b ^ k^^). It can be shown that in the case on 
an inviscid liquid the wave resistance per unit length corresponding to the 
external pressure distribution (g) is given by 



R 



Pr 



kie-^^^' + k2e-^^^A {V > c™") (3) 



(remember that for an inviscid liquid, R = for V < c*"*"^). In Eq. (|]), the 
wave numbers ki and ^2 denotes the two (real) solutions of {g/k + 'yk/ pY^'^ = 
V (see Eq. (|1])). A brief inspection of Eq. (^ shows that the wave resistance is 
a decreasing function of the perturbation velocity V. In the limit V ^ c'"*"', 
Eq. (|) reduces to i? ~ {Pq /'y)e~'^^^^'^~ )(V'/c™»") ^ ^.j^g velocity V decreases 
towards c™'", the wave resistance Eq. (^ becomes unbounded. This result is 
directly related to the fact that as V approaches c™", the energy transferred 
by the moving pressure distribution cannot be radiated away. 

We now turn our attention to the case of a viscous liquid and investigate 
how the wave resistance (0) is modified by the liquid viscosity. In order to 
calculate R, we may imagine a rigid cover fitting the surface everywhere, as 
suggested by Havelock The assigned pressure system P^xt is apphed to 
the liquid surface by means of this cover; hence, the wave resistance is simply 
the total resolved pressure per unit length in the x-direction |^ : 

R = -J Pextix) (^C(a;)^ dx (4) 

where ({x) denotes (in the frame of the perturbation) the displacement of the 
free surface from its equilibrium position. Let ( and Pgxt denote the Fourier 
transforms of ( and Pext, respectively. Using the Navier-Stokes equation for 
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a viscous fluid along with the appropriate stress condition at the free surface 
[^], one can relate ( to Pext through the relation 



(2z/A;2 - iVkf + g \k\ +^\kf - 4u^\kf]jk^ - i^k 



C = -SPext (5) 
P 



where the parameter v is the kinematic viscosity of the liquid. Inserting 
Eq. (H) into Eq. @ we obtain 

P2 / 7 ik^ exp (-A^v^ k) 

7^7 \{ (2eo vk^- iky + \v-*k + k^- Ael k^^k^ - 

(6) 

where v = V/c"^^^ and eo = t^c"^*"'/^ . The symbol 3fJ represent the real part 
of a complex expression. The wave resistance Eq. (P) is shown graphically 
in figure (1) as a function of the reduced velocity v for cq = 3 10^^ and 
b = 2.5 10~^ K,~^ (the inset highlights the behaviour of R at low velocities). 
Two important features should be noted in comparison to the inviscid case: 
(a) First, while R increases steeply near V = c™", it remains bounded; (b) 
Secondly, as soon as the perturbation velocity V is greater than zero, the 
wave resistance takes finite values. This late result is a direct consequence 
of the internal viscous dissipation inside the liquid. 

For eo much smaller than unity, Eq. (P) can be simplified and the above 
two features can be recovered analytically (for water with 7 = 73 mM m~^ 
and p = 10^ kg m~^, eo ~ 3 10~^ <^ 1). Using standard mathematical 



technics ||10[, one can show that the wave resistance displays a maximum of 



p2 

R Rmax ~ , — ("7) 



7 V^o 
for 

V = V^ax ~ C"™ (1 + eo) (8) 

On the other hand, in the limit V -C c™*", the wave resistance R varies 
linearly with the perturbation velocity: 
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This linear behavior can be observed in the insert in figure (1). 

Let us conclude by a few remarks. In the calculations made above for 
the wave resistance R, we have used Rayleigh's linearized theory of capillary- 
gravity waves 0. We have shown that one of the effects of viscosity was to 
cutoff the unbounded response of the liquid predicted by the inviscid model 
as \ c"^™. For a small viscosity, the response of the liquid remains, how- 
ever, rather large near c™" (see figure (1)), and it would be of some interest 
to take nonlinear effects into account in the calculation of R. This is beyond 
the scope of the present letter. For a recent review of nonlinear capillary- 
gravity waves, the reader is referred to the work of Dias and Kharif |Tl|]. In 
the present study we have emphasized the asymptotic behavior of the wave 
resistance for a liquid of low viscosity. We hope to explore the high viscosity 
limit in a subsequent report. Note also that the calculations presented in 
this letter assumed an external pressure distribution localized along a band. 
Further work will assess the effect of viscosity for a pressure field symmetrical 
around a point. 
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Figure caption : R* = tt'j Pq^R as a function of f = V^/c™", with R 
being the wave resistance Eq. (P), eo = 3 10"'^ and b = 2.5 10~^ k,~^. The 
inset shows the behavior at low velocities. 
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